	PRO helicity_globals,dir
; Usar leer_globals.pro

leer_globals,dir,g,NN
t=fltarr(NN+1)
hel_p=t ; proton energy
hel_e=t ; magnetic energy
hel_c=t ; electron energy

t(*)=g(0,0:NN)
hel_p(*)=g(12,0:NN) ; proton helicity
hel_e(*)=g(11,0:NN) ; electron helicity
hel_c(*)=g(7,0:NN) ; cross helicity (<u_perp.b_perp>)

!p.multi = [0,1,2]
max_y = 1.2*max([max(hel_p), max(hel_e)])
min_y = 1.2*min([min(hel_p), min(hel_e)])

window,0,xs=800,ys=800
plot,t,hel_p,linestyle=0,xtitle='Time',tit='Proton Helicity',xr=[0,max(t)],yr=[min_y,max_y],xsty=1,ysty=1,charsize=2.0
plot,t,hel_e,linestyle=0,xtitle='Time',tit='Electron Helicity',xr=[0,max(t)],yr=[min_y,max_y],xsty=1,ysty=1,charsize=2.0

!p.multi = 0
window,2,xs=650,ys=650
plot,t,hel_c,linestyle=0,xtitle='Time',tit='Cross Helicity (<u_perp.b_perp>)',xr=[0,max(t)],yr=[1.1*min(hel_c),1.1*max(hel_c)],xsty=1,ysty=1,charsize=2.0
end
